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■ We study the scale at which one can generate the lepton asymmetry 

of the universe which could then get converted to a baryon asymmetry 
during the electroweak phase transition. We consider the possibility 
that the Yukawa couplings are small but sufficiently large to generate 
enough lepton asymmetry. This forbids the possibility of the (B — L) 
breaking scale being the electroweak scale. 
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1 Introduction 



In most grand unified theories (GUTs) the baryon asymmetry of the universe 
is generated during the GUT phase transition [1], 0, §. In these models the 
generated asymmetry also implies an equal amount of lepton asymmetry 
and hence there is no net (B — L) asymmetry. On the other hand if the 
electroweak phase transition is a second order phase transition, then any 
primordial (B + L) asymmetry generated during the GUT phase transition 
will be washed out [4| . 

This situation can be saved in models where (B — L) is broken at some 
intermediate scales. In this case a (B — L) asymmetry can be generated 
through higgs decay or heavy Majorana neutrino decay if there is appropri- 
ate CP violation |J-[11]. The out-of-equilibrium condition then imposes a 



lower bound on this symmetry breaking scale to be around 10 7 GeV. This 
bound is dependent on the fact that the Yukawa couplings are larger than 
10~ 5 . Although esthetically this number sounds reasonable, nothing tells us 
definitely that the Yukawa couplings cannot be smaller than this. For exam- 
ple, if the Yukawa couplings relating the left handed leptons to the first gen- 
eration right handed heavy neutrinos are similar to the Yukawa couplings for 
the electron, then the out-of-equilibrium condition can be satisfied for even a 
TeV scale for left-right symmetry breaking. But the same Yukawa couplings 
enter in the expression for the generated {B — L) asymmetry, which may 
then be very small. In this article we study systematically the Boltzmann 
equations for the generation of the lepton asymmetry and hence (B — L) 
asymmetry to find out the lowest possible left-right symmetry breaking scale 
which satisfies the out-of-equilibrium condition and generates enough baryon 
asymmetry after the electroweak phase transition. 

In the next section we briefly review the leptogenesis scenario, where one 
generates a lepton asymmetry when the right handed Majorana neutrinos de- 
cay. This then gets converted to a baryon asymmetry during the electroweak 
phase transition. Subsequently, we discuss the Boltzmann equations and the 
possible solutions for a low energy left-right symmetry breaking. 
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2 Model for leptogenesis 



It was first proposed by Fukugita and Yanagida || that in extensions of 
the standard model, which include singlet heavy right handed neutrinos, 
it is possible to generate a lepton asymmetry at some intermediate scale 
(~ 10 10 GeV), which can then get converted to a baryon asymmetry during 
the electroweak phase transition. The source of CP violation is in the mass 
matrix of the neutrinos. The out-of-equilibrium condition is satisfied with 
small Yukawa couplings. In this scenario the heavy right handed neutrinos 
decay to light left handed neutrinos in out-of-equilibrium. The amount of 
asymmetry thus generated depends on the amount of CP violation. 

In this approach it is expected that at the time of the decay, there are 
enough right handed neutrinos given by their equilibrium distribution. How- 
ever, in practice it is difficult to achieve this distribution ||. Even though 
we assume that gravitational fields have generated equilibrium numbers of 
right handed neutrinos, during the period 10 18 GeV — > 10 10 GeV the equilib- 
rium distribution cannot be maintained since the right handed neutrinos are 
sterile under the standard model and as the universe expands the number 
density reduces. When these fewer right handed neutrinos decay in out-of- 
equilibrium, they cannot generate enough lepton asymmetry. 

It has been pointed out that if these right handed neutrinos take part in 
any gauge interaction, then the decays of the gauge particles can produce an 
equilibrium distribution of the right handed neutrinos. The left-right sym- 
metric models are the natural extensions for this, where the heavy neutrinos 
interact with the right handed gauge bosons strongly ||. As a result, the 
decay of the right handed gauge bosons can generate an equilibrium distri- 
bution of heavy neutrinos. In addition to the heavy neutrino decays, lepton 
number violating higgs decays also contribute to the generation of lepton 
asymmetry in these models M. 

In the following we shall thus consider a left-right symmetric extension 
of the standard model. We consider the symmetry breaking chain, SU(3) C x 
SU(2) L x SU(2) R x f/(l) (B -L) [= G LR \ ^ SU(3) C x SU(2) L x U(1) Y [= G M ] 
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w SU(3) C x U(l) em . The symmetry breaking Glr — > G st d takes place when 
the right handed triplet higgs field Ar = (1,1,3,-2) acquires a vacuum expec- 
tation value (vev). In this model (B — L) is a local symmetry. The breaking 
of the group Gl/j also implies spontaneous breaking of (B — L). Left-right 
parity implies the existence of another higgs field Al which transforms as 
(1,3,1,-2) under G LR . A higgs bi-doublet field (1,2,2,0) breaks the elec- 
troweak symmetry and gives masses to the fermions. 

The fermion content of the model is, q iL = [3, 2, 1, 1/3], q iR = [3, 1, 2, 1/3], 
£iL = [1,2,1,1] and £r = [1,1,2,1], where % = 1,2,3 corresponds to three 
generations. The right handed neutrinos (iVj = vm) are contained in im 
and we do not have to include them by hand. The Yukawa couplings in the 
leptonic sector are given by, 

£y«fc = fij^iL^jR<t> + fLijtiL C ijl J A\ + fRij£iR C £jRA R . (1) 

The scalar potential has many more terms compared to the standard 
model. We write down only those terms which contribute to the generation 
of the lepton asymmetry of the universe, 

C int = giAlA^ + A L A^00) + h.c. (2) 

The vevs of these fields are not independent. We consider the minimum of 
the complete potential which satisfies vl vr and vl ~ v 2 /v R , where v^r 
and v are the vevs of the fields Al,r and <fi respectively. We also assume 
that the left-right parity (D— parity) is not broken, and hence the masses of 
the fields and Ar remains the same even after the breaking of Glr, i-e., 
m A L = m A R = m A ~ Vr. At this scale vr, the (B — L) local symmetry is also 
broken by two units, which gives rise to Majorana masses of the neutrinos 
and neutron-antineutron oscillations. 

The A L R can now decay into two leptons, while A L R decay into two 
antileptons: 

A L , R -> £ l ,r + £l,r- (3) 

AU - ti,R + tl,R- (4) 
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These interactions, along with the scalar interactions 

- f + f , 

give rise to lepton number violation. The interference of the tree level di- 
agram and the one loop diagram of Fig. 1 can then give rise to a lepton 
asymmetry in the decay modes of @ and @ for the left-handed triplets 
given by, 

where F(q) = ln(l + l/q 2 ). The quantity [g*flijfikfjk] contains a CP violat- 
ing phase and so can be complex. 

Although the left handed higgs triplets Al can generate lepton asymme- 
try through their lepton number violating interactions, the scattering process 
Al + Al — > Wl + Wl will make their number density the same as the equi- 
librium density below the GUT scale so, at low energies, they cannot satisfy 
the out-of-equilibrium condition. So for the generation of lepton asymmetry 
at low energy the main contribution comes from the heavy right handed neu- 
trino decays. The vev of A# spontaneously gives a Majorana mass to the 
right handed neutrinos. This in turn allows the decay of Vr into a lepton 
and an antilepton, 

n % -> e jL + $, (6) 

- £ jL c + <f>. (7) 

In the case of decays of the right-handed neutrinos there are two types of 
loop diagrams which can interfere with the tree level decays of (^) and (0) 
which are shown in Fig. 2. The interference of the tree level diagram and 
the one loop diagrams of Fig. 2 generates a lepton asymmetry given by, 

^^T-^MMilfJJJi}^- (8) 

^\jik\ JRkk 

In addition to these one loop vertex type corrections, there are self energy 
type corrections [UJ which also contribute to the lepton asymmetry of 
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the universe (Fig. 3). This contribution can also be viewed as a new indirect 



CP violation entering in the mass matrix of the heavy neutrinos [|TTJ. This 
contribution becomes much stronger than the other contribution when the 
heavy neutrinos are almost degenerate. 

When the mass squared difference between the two generations of heavy 
neutrinos are large, the interference of the tree level and the one loop di- 
agrams of Fig. 3 has already been calculated in the literature ||10|| . This 
contribution has also been calculated in the other approach, where this is 
viewed as a new indirect CP— violation 5, which enters through the mass 
matrix. Both the results come out to be the same and that confirms the 
equivalence between the two approaches. However, the latter approach has 
the advantage that it allows us to calculate the amount of CP violation in 
the limit when the mass squared difference is very small. 



In the case of a small mass difference, S reads |L1 



6 = 1*4*0 f M \ (9) 
y Ml -Ml 

where 



C = -ilmEOSi/oa) E fpiff*)\ L I 12 + ^ i 1 , ia ] ( 10 ) 



77 a p \Sa I /oil 2 ^2a \f< 



a2\ 



This contribution becomes significant when the two mass eigenvalues are 
close to each other. It indicates a resonance like behaviour of the asymmetry 
if the two mass eigenvalues are nearly degenerate. For very large values of 
the mass difference the two contributions e v and 5 are of the same order of 
magnitude. 

We now have to consider only the Al, Ar and Ni^ decay processes. We 
assume N% to be the lightest of the right handed neutrinos. If the masses of 
N\ and N2 are almost degenerate, their decay widths can become larger than 
the mass difference. In this case both the neutrino decays will contribute to 
the lepton asymmetry of the universe. The decay widths for Al ,r and Ni 
are, 

Ta lr = 1 -^^M a and Y Ni = ^M Nl (11) 

lD7T l07T 
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where Mjy is the mass of N±. Since we assumed Mn < Ma, at low energy 
the Al decay will erase all lepton asymmetry and then the Ni decay will 
generate the required asymmetry. For this reason while working the details of 
the Boltzmann equation we take the effect of only e = e u and 5. We shall now 
proceed to solve the Boltzmann equations including all these contributions 
and the scattering processes. 



3 Solutions of the Boltzmann equations 

The evolution of lepton and neutrino densities is governed by the Boltzmann 
equations. We start by deriving the Boltzmann equation for the neutrino 
number density n^. The equation governing the evolution of rii is:0 

x {-f*(Pi)\Mo\ 2 + \{l + e)|7W | 2 //(Pi)/^(P2) 
+ ^(l-e)\M \ 2 f I (p 1 )U(p 2 : 

= |dn 4 dn 1 dn 2 (27r)V 4 )(j9 t -Pi-p 2 ) 

{-fM + f?(Pi)} IM.I 2 + 0(e, fi/T) 
= -Tiim-n?), (12) 

where n^ q is the equilibrium number density of the Ni, and Tj is the thermally 
averaged decay width of iVj. The term on the left-hand side is the time 
derivative of rij, plus a term which accounts for the dilution effect of the 
expansion of the universe. The integration is over phase space dll and the 
phase space densities /, are given by the Maxwell-Boltzmann statistics: 

exp[-^^], f t (E) = exp[-^4], U{E) = exp[-^±^], 

(13) 

where [i is the chemical potential. The matrix element .Mo is defined by: 
\M(N^U)\ 2 = \M(l$ c ^N)\ 2 = ^(l+e)\M \ 2 , 
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\M(N ^ 1<5> C )\ 2 = \M(l$ -> N)\ 2 = -(l-e)\M \ 2 , (14) 

where e is a measure of CP- violation. 

It is more convenient to work with the variables: 

Yi = m/s, x = Mi/T=[2H(x = l)t]*, (15) 

where Mj is the mass of Nf, s = g*n^ is the entropy density of the universe; 
g* is the total spin degrees of freedom; n 7 is the equilibrium photon density of 
the universe and Y$ the number of neutrinos per co-moving volume element. 

Thus: 

^ = -K 1 x(Y l -Yn, (16) 



where we have used: 



and: 



n 7 = s/g*, ^ = ~ 3s ^ = SsH, (17) 



K = ^ = 1} 7= ^ (18) 

H(x = iy 7 r,(x = i)' 



with: 



veq _ eq/ _ J a: 1 X<1 



In solving the Boltzmann equations, we make the further change of variables, 
X = g*Y, thus: 

( ^ = -K 1 x{X l -XD. (19) 

This is the Boltzmann equation for the evolution of neutrino number density, 
with the initial condition JQ(0) = I. 

We now derive the Boltzmann equation for L — |(Z — I), where we have 
to take into account the processes I + $ c <-> 1+ $ mediated by a right-handed 
Majorana neutrino, as well as the processes: 

l^N$ c . (20) 
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The Boltzmann equation for the number density of the light left-handed 
leptons is: 



The origin of the various terms is described below. 

The first interaction term comes from direct CP- violation in the decay 
of the heavy Majorana neutrinos and is proportional to e, the CP-violating 
phase and the squares of the matrix element \M \- The second term describes 
the generation of leptons from an initial neutrino state \Ni >, with number 
density rii. This converts into an antineutrino of a different generation \N? > 
which then decays into a lepton and a higgs <£> with decay width: 



The last term takes into account the 2 <-> 2 scattering mentioned above. The 
corresponding equation for I is obtained as usual, by interchanging I <-> I, 
e <-> — e, etc. To obtain the Boltzmann equation for the evolution of rii = 
\{n\ — nj) we subtract the equation for nj from that for rii and multiply by 
a factor of 1/2: 

h L + 3Hn L = eYiirii — n^) — n L {n e ^ ' jn^)Y i — 2n L n y < a\v\ >, (23) 

where \M'(l$ c — > I$)| 2 and \M.'(l$> — > /$ c )| 2 are the squares of the matrix 
elements for 2^2 L-nonconserving scatterings with the part due to real, 
intermediate-state TV's removed. Here the quantity: 



< a\v\ > = j dn 1 dn 2 dn 3 dn 4 (27r) 4 (5 {4) (p 1 +p 2 -P3 -Pa) 

fi(pi)fi( P 2)\M'(l$ c ^m\ 2 K (24) 



■hi + 3Hn t 



Jdu N dn x dn 2 (2tt) 4 5^ ( PN - Pl - P2 ) 

x [-(1 - e)f l i Pl )Uc(p 2 ) + (1 + e)f N (p N )} \M \ 2 
+2 JdYli dU 2 dn 3 dn 4 (2tt) 4 5 w ( Pl +P2-P3- Pa) 
x -fi(pi)Uc(p 2 )\M'(l$ c ^l$)\ 2 
+f T (pz)U(pA)\MXU^l<$> c )\ 2 ] + 0(en l + n 2 ). (21) 
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is the velocity-averaged 2 <-> 2 L-violating cross-section. The presence of the 
term — er^ 9 is due to the CP-violating part of \M'(l$ c — > /$)| 2 - \M'(l$ —> 
/$ c )| 2 . 

In parallel to the calculation for the n^, we obtain: 

= ei^^z - K") - 3,^*7* - , (25) 

where T s = n 7 < cr|u| >, with the initial condition Yl(0) = 0. 

To solve the set of coupled differential equations, we take for I\, 7 and 

r s : 

III • ^ 

and: 

r, = ^f!, *<J . (28) 

The functions 7, T,, and T s must be specified in the region x ~ 1, in order 
to find a solution for the differential equations. We take: 

7 = 1 -exp(-x), (29) 



and: 



167T 

3/ 4 M, 



1 _ exp(-x)) , (30) 



(l-exp[-2/(3* 2 )]), (31) 



Attx 

which are a good approximation to the functions fl2"6|), Q2"7| ) and ( J28|) above, 
in the regions iCl and s>l. The numerical solutions to the Boltzmann 
equations are then obtained for different values of /, K and e. We summarize 
our observations below. 

For K <C 1, the amount of lepton asymmetry grows cubically to a con- 
stant asymptotic value we call it Y^ sym (K <C 1) as shown in figure 4. In all 
these figures we have taken some representative values for the couplings to 
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be of the order of 1. The nature of the curve is independent of the choice of 
the Yukawa couplings. This asymptotic value is given by, Y^asym) = 
This is the case when the decay rate of N\ is less than the expansion rate 
of the universe. In this case the scattering rates are also less than the ex- 
pansion rate of the universe. If the mass difference between N\ and N 2 is of 
the same order of magnitude to their masses, the contributions of e and 5 
become comparable. This case has been discussed extensively in the litera- 
ture and the constraint on the scale of (B — L) breaking obtained from this 
condition is Mr > 10 7 GeV with the Yukawa couplings to be of the order of 
10~ 5 . Somewhat smaller values of the Yukawa couplings can reduce the scale 
of (B — L) breaking to a lower value. However, then the amount of lepton 
asymmetry will also be inadequate, unless there is large hierarchy among the 
Yukawa couplings of different generations. 

When the mass difference \M 1 - M 2 \ ~ 10~ 3 Mi j2 , 6 can be three to four 
orders of magnitude larger than e. In this case we can consider |/| ~ 10~ 7 
and still get adequate amount of lepton asymmetry. This will then allow us 
somewhat smaller right handed symmetry breaking scale M# > 10 5 GeV. 

When K is of the order of unity or more, the lepton asymmetry vanishes 
before T = M\. From T = M\ onwards the lepton asymmetry starts in- 
creasing from its initial value of Yl = 0, but as it approaches the asymptotic 
value of Y^ sym (K -C 1) the scattering processes becomes important and start 
depleting it exponentially (as shown in figures 5 and 6). Unlike the common 
folklore that when the system is in equilibrium the asymmetry falls exponen- 
tially fast, soon the scattering processes become unimportant and the lepton 
asymmetry reaches its new asymptotic value, which is less than — . For 
K = 1 the suppression factor is about 1/5 and for K = 5 it is .02. Figure 7 
shows the fall of this asymptotic values of the lepton asymmetry for different 
values of K. For K = 1000 the suppression factor is about 8 x 10 -6 . In grand 
unified theories where the heavy gauge bosons decay generates an (B + L) 
asymmetry it was shown that this suppression factor is ~ [^(lnii') ' 6 ] -1 . 
But in the case of leptogenesis, as figure 7 shows, this suppression is much 
faster than linear. For large mass difference \M\ — M 2 |, Yl is approximately 
proportional to f 2 when all the /'s are of the similar order of magnitude 
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and K ~ lO -3 / 2 ^. Since K and Yl both are proportional to \ f\ 2 , we can- 
not improve the bound on the right handed symmetry breaking scale with 
large K. However, since the suppression for large K is almost quadratic we 
can at most lower the scale of left-right symmetry breaking by one order of 
magnitude to about 10 4 GeV. 

For very large K, if \f\ 4 > 10 _17 M 1 , the scattering processes become 
larger than the expansion rate of the universe. In this case the lepton asym- 
metry decreases exponentially and never reaches any asymptotic value. In 
this region the actual equilibrium condition is satisfied. 

4 Summary and Conclusion 

We have studied the possible scale of left-right symmetry breaking, which 
can generate enough lepton asymmetry of the universe. For simplicity we 
assumed that the order of magnitude of the Yukawa couplings of the heavy 
neutrinos are similar. It was shown that the right handed symmetry breaking 
scale can be lowered to 10 4 GeV at most. Slight deviation from the out-of- 
equilibrium condition can still generate similar amount of leptons asymmetry, 
but that cannot change the scale of the symmetry breaking. 
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Figure 1: Tree and one loop diagrams of lepton number violating triplet higgs 
A L decay 
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Figure 2: Tree and vertex correction type one loop diagrams contributing to 
the generation of lepton asymmetry 
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Figure 3: Tree and self energy correction type one loop diagrams contributing 
to the generation of lepton asymmetry 
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Figure 5: For K — 1, lepton asymmetry starts depleting before reaching a 
constant value. The asymptotic constant value is thus much less than e. 
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Figure 7: The asymptotic value of the lepton asymmetry for different values 
of K in log w — log w graph. For K = 1000 the lepton asymmetry drops to 
8 x 10~ 6 . 



20 



